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^ ABSTRACT 



' We consider a scalar massless quantum field model, at finite temperature T, both renormalizable and asymp- 
, totically free. Focussing on the singular structure of the effective perturbation theory about the light cone, 
C3 ' several new insights are put forth, regarding the interplay between hard thermal loop resummation and the 
overall compensation of collinear singularities. 
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1. Introduction 

The inherent non perturbative character of finite temperature quantum field theories has been recognized 
a few years ago on the basis of general group theoretical arguments [1]. At the same time, the formal 
perturbative series themselves were shown to necessitate a full redefinition of their original form so as to 
yield (hopefully) sensical results. This is the resummation program devised by Braaten and Pisarski and 
also by Frenkel and Taylor [2], hereafter referred to as "effective theory" for short. 

In order not to be just an academic recipe, such a program implicitly requires that the original and/or 
re-organised perturbative series be infrared safe. In the recent years, a great deal of efforts has been devoted 
to the study of the infrared structure of the original perturbative series, with the conclusion that, roughly 
speaking, the situation at non zero temperature (T ^ 0) was not worse, globally, than at zero temperature 
(T = 0)[3]. 

But re-defining a perturbative series is rarely trivial an operation. Properties which were known to hold 
true for the original series may become much harder to control in the redefined one. In our opinion, the 
collincar singularity recently found in thermal QCD by using the effective perturbative series [4,5] might 
reveal a salient illustration of this fact. Note that we are aware that this point of view differs, at first sight, 
from the ones adopted in the solutions proposed in [6] and [7]. Very recently, though, it has been realized 
[18], that some unexpected connection could exist with the analysis proposed in [6]. In this paper we will 
restrict ourselves to a scalar model and analyse the singular structure of the effective theory in much details. 

A preliminary report of the present work appeared in [8], which is here corrected, completed and enlarged. 
The paper is organized as follows. Section 2 is an introduction to the model, and the necessary elements 
of the particular real time formalism which is used, are given. In section 3, the singularity structure of 
a whole series of diagrams is investigated, deferring some lengthy details to an appendix. The results of 
this section rely on an approximation which consists in keeping only the potentially most singular part of 
the hard thermal loop (HTL) self energy. This approximation gets completed in section 4, where the other 
terms are taken into account. Further consequences are drawn in sections 5 and 6, respectively related to the 
problem of the T = and T ^ singularity mixing, and to the possibility of an unambiguous renormalisation 
constant definition in a thermal context. Eventually, a conclusion gathers our results. 



2. The model 

In order to get rid of unessential complications, we deal with the toy model provided by a massless hermi- 
tian scalar quantum field. It is endowed with a cubic self interaction, of coupling strength g, rcnormalisable 
at D — 6 space-time dimensions. Furthermore, its quanta are assumed to form a plasma in thermodynamical 
equilibrium at temperature T. Infrared as well as ultraviolet singularities are taken care of by working at 
D = 6 ± 2e space-time dimensions respectively. Also, we will use the convention of upper case letters for 
D-momenta and lower case ones for their components, i.e, P = (po,p). 

Use will be made of the material and results of a previous work [9] where a calculation performed at second 
non trivial order of bare perturbation theory, 0(g 4 ), was able to display the overall compensation of infrared 
and collinear singularities for some particular topology. 

The process under consideration consists in a "Higgs" particle of D-momentum Q — (q, 0) in the rest-frame 
of a plasma in equilibrium at temperature T. But, contrary to the case studied in [9], this Z?-momentum is 
here assumed to be soft, that is q = O(gT). Also, our analysis will take place in the framework of the real 
time formalism, within "the basis" provided by the retarded and advanced free field functions, Ar(K) and 
Aa(K) respectively [5]. At any step, though, we have checked that the same expressions are obtained in the 
more customary real time formalism of the Kobes and Semenoff rules (the "1/2" formalism, say). 

The real part of the thermal one loop self energy reads [9] 

ReS T (P) = - ^ / dfc k D - 3 n(k) V P / dx , 1 , , ; (2.1) 

rv ' (4^)3+-r(2 + e) Jo J-i P 2 + ZHwo ~ px) K ' 
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where a combinational weight factor 1/2 is taken into account, and where P stands for a principal value 
prescription. Also, throughout the calculations, use will be made of the following relations for statistical 
factors 

n(k Q ) = (e^" - I)" 1 , n B (k Q ) = (e^ - l)" 1 , e(k a )(l + 2n(fc )) = 1 + 2n B (fc ) (2.2) 

The phase-space domain where collinear singularities come from has been recognized to be determined by 
the condition P 2 « p 2 « T 2 . In this domain, the function {PeEy}, that we choose to describe in terms 
of the two independent variables P 2 and p 2 , is approximated by a real function A, defined for all P 2 values 

A(P 2 y) = -m 2 \ J|l - i±|^ \^\ (0(P 2 ) + cos (7 r £ )e(-P 2 ))| (2.3) 

where to 2 , the thermal mass squared, is given by 

™ 2 = 2 (£^3 C(2 + 2e)T 2+2£ (2.4) 

Expression (2.3) displays the full HTL (that is, leading, order g 2 T 2 part) entailed in the real part of the 
thermal self energy function. Only for space- like values of its argument P 2 , does the thermal self energy 
function develop an imaginary part of same leading order g 2 T 2 . This corresponds to the so called "Landau 
damping mechanism" , and one can write eventually 



v r = A+iu . b(7>V) - ----- ^r^i-p 2 ) 1 - 1 ^ ^ (->■■» 



where Tip denotes the Feynman self energy. Indeed, the functions A and B can be seen as obtained through 
the usual Feynman prescription {pq — ► po + iwpo}, out of one analytical self energy function £ defined in 
the complex energy plane 

A = 1im - (S(p + iVPo,p) + ^(Po - iVPo,p)) 
B = lim 7T- ( s (Po + iVPo,p) - S(p - iVPo,p)) 

r]=0 2,1 

with E the light cone approximate self energy 



1 P 2 l + 2 1+e e (P 2 

1 ~p^\ X 4^ 




^ = -m 2 --T{^-^i^[-) \ (2-6) 



Note that in (2.6), the limit e = can be safely taken and that the thermal self energy function is a regular 
function of P 2 in any neighbourhood of the light cone. One gets effectively 

S(P 2 ,p 2 ) = -2™ 2 ^filn(^)-l") (2.7) 



p 2 \2 y P 2 

which, expressed in terms of the variables P 2 and p 2 , is nothing but the more familiar form [12] of the HTL 
self energy of the model, with its second kind Lcgcndre function Qi{po,p) approximate form at \po\ ~ p. 
In the R/A real time formalism we will be using, retarded and advanced functions must be used. They are 
simply related to the Feynman ones by the set of relations 

Sii(p) = e(p )s F (P) + e(-po)EMP), e a(p) = e(po)sj.(P) + e(- Po )z F (P) (2.8) 

In the region V where one has P 2 « p 2 « T 2 and 0{p^) = gT, effective propagators *A Q (P) must be 
used so as to get the full leading order g 2 T 2 correction to the free quantities, with 

• A °< f > = P*-s„ (P W —*•*■ (2 - 9) 
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A point to be stressed is that the domain which requires that effective perturbation theory be used, is also 
the domain where mass (or collinear) singularities stem from, thereby opening up the very "window" where 
the light cone singular behaviour of the effective thermal expansion must be studied. 

An important remark is in order. Admittedly, there is no HTL in the three-point function of the (g<p 3 )6 
theory [10], in the sense that g 2 T 2 leading parts, eventually, cancel out, leaving order g 2 T 3 / 2 pieces at 
most. However, in [10], not all the terms were considered and, more generally, it should be noted that this 
statement is not sufficient for ignoring the one loop three point function. For example, it has been stressed 
recently [7] that, over collinear configuration of external soft momenta, a special enhancement of bare orders 
of magnitude, results out of collinear singularities, leading to a breakdown of the HTL rcsummation program. 
However, besides the new internal difficulties encountered by the improved HTL resummation program itself 
[24], it is important to realize that such a situation docs not show up here. This is basically because the Higgs 
particle has timelike D-momentum Q = (q, 0). Explicit calculation then shows that no collinear singularity 
ever develop in that case. Now, a wider analysis of the three point function is certainly worth undertaking, 
but falls beyond the scope of the present article where one loop vertex corrections will not be considered. 



3- Singular structure 

We begin with considering the case of N one-loop self energy insertions along the P-line as depicted 
on Fig.l, with N' = 0. Then, the retarded "polarisation tensor" n^2(Q) of the" Higgs particle" (hereafter 
written TI^\q) for short) can be shown to admit the following expression (see (61) in [5]) 

nf(Q) = -i / ^(l + 2n(p ))msc P (A { R N \p)A R (P')V 2 RA (P,Q,-P')) (3-1) 

where Vrra is the bare vertex with two external lines of the retarded type, corresponding to D-momenta 
P and Q, and one advanced external line of D- momentum —P'. The convention for the flow of momenta is 
that their sum vanishes: P + Q + (-P') = 0. At zcroth order, Vrra is simply given by the "electric charge", 
e, which couples the Higgs particle to the hermitian scalar field. Eventually, the prescription "Discp" means 
that the discontinuity in the variable po is to be taken. In the {R/A} real time formalism we are using, this 
is simply achieved by writing 

Discp Frps(P, Q, R) = F Rps {P, Q, R) - Fa(38{P, Q, R) (3.2) 
with (3 and S any of the two R, A possibilities. 

In above equation (3.1), denotes a N one-loop self energy corrected free field function, and n^((Q), 

the ensuing corrected tensor. Omitting unessential factors, its imaginary part reads 

ImnW(Q) = J d D P(l + 2n( Po )) £ (^) ( 5(P' 2 ){^ TC ( Po )^ jV) (P 2 )Rc(Sg(P)) - P ^fffjf^ > ( 3 - 3 ) 

where e(p ) is the sign of p , and the N th derivative of the ^-distribution. In view of equations (2.3), (2.5) 
and of their bounded P-domain, note that the real and imaginary parts of the {£|> TL (P)} ^-insertions are 
acceptable test-functions, yielding well defined expressions when acted upon with 5^ N \P 2 ) and P / (P 2 ) N+1 
distributions respectively. How things are working out is instructive to observe. At N = 1, we have from 
the (^(P^-term 

— e 2 f f dP 2 

WY^J d ^ j /^(l + 2n(p o) ) e ( Po ) 

x £ (p' )S(P' 2 )7rS(P 2 )-^A(P 2 ,p 2 ) (3.4) 
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where \po{P 2 ,P 2 )\ = \J P 2 + P 2 ■ Substituing (2.3), we get immediately 

where fie is the total solid angle of the model. The singular piece coming from the term involving the 
principal value, more involved, is 

D-l 



1+2l+££ ^/^(l + 2n( P )) £W (0 



(2tt) d - 1 4 e 

dP 2 / 2 P 2 



where (2.5) has been used. By integrating over P 2 first, which yields a £ 1 mass singularity, and then over 
p, one gets 

I ne2m2 + 2n B f^W-) 1 ± 2l+££ Sin(7r£) f 3 7) 

The sum of the two terms displays the compensation of the collinear singularity, as noticed long ago by many 
authors for scalar as well as for gauge theories. 

At N = 2, one has ReT, 2 R (P) = A 2 {P) - B 2 {P). The function B is of the order of (-P 2 ) 1+e with e positive, 
and therefore does not contribute when acted upon with the 5^ 2 \P 2 ) distribution. Indeed, for exactly the 
same reason, only the {^^j-part of ReS^(P) will ever contribute at general N, because of the 8^ N \P 2 ) 
distribution. In the present case, N = 2, one gets 

Sin 2 / _ n ,9,\ ,m 2 



(' +-*<!>) <t> 2 <«> 



(2TT) - 1 Aq q 

which, by recalling that q is of the order of gT, is consistently seen to be on the same order as (3.5) and 
(3.7) at N = 1. Now, the part involving the principal value distribution reads 

e 2 1 + 2 1+e e sin(Tre) 



/ d °~ lp I ^( 1 + Mpo)Mp' )5(p' 2 ) 



{2tt) d - 1 ¥ ire J " J 2 Po v v™//~v™/~v > y (p 2 y 

2n2 i j ol+e,- o^/V,-^ E>2 



„. . / m 2 P 2 , l + 2 1+£ £ simVeVP 2 , \ / ,. P\ 1+F \ ^, ^ ,„ ^ 
x (1 - — p ^1^1 £ cos( 7 r £ ))J (-™»(_) 1+e J 9(-P 2 ) (3.9) 



Integrating out over P 2 and p, one gets 



e 2 D 2/ „ R/ q,\ l + 2 1+£ £ sin tte ). m 4 ./l 1 + 2 1+£ £ , ,\ . 

+71" 7 — 7n— r — - 1 + 2n s (- i — -( cos(tt£) 3.10 

(2tt) d - 1 Aq q V V 2V 4 £ tte V £ ' \s A e 2e y ') V ' 

Expanding in powers of e the terms 

l + 2 1+£ £ sin(7T£) (\ l + 2 1+£ £ , A 1 . 2xN 

— ^ — 7^ cos(vr£) = — (1 + 0(e 2 )) 3.11) 

4 £ 7T£ \£ 4 e 2£ V 7 2£ V V ^ V ^ 

the sum of the two terms (3.9) and (3.10) is readily seen to display the expected singularity compensation. 
Indeed, inspection of the first values of N (1 < N < 5), shows that one may drop the terms 2 £ , 4~ e appearing 
in (2.3) and (2.5) without spoiling at all the singularity cancellation mechanism. In what follows, we will 
take advantage of that simplification in order to avoid overcharging with details. 

At N = 3, one has 

T? R (P 2 ,p 2 ) = (A 3 - 3AP 2 )(P 2 ,p 2 ) + ie( Po )(3A 2 B - B 3 )(P 2 ,p 2 ) (3.12) 
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Of the real part, only the {A 3 } term will contribute according to the remark above. The resulting singularity, 
the stronger, of order £~ 3 is given by 



((-l) 3 /3!) ne( Po )S^(P 2 ) A a (P 2 ,p 2 ) = ns(p )5(P 2 ) {-m 2 /ep 2 f (3.13) 

where, for short, we do not write the remaining integration over p (see the remark after (3.20)). The 
imaginary term {— B 3 } is easily integrated over P 2 and yields 

, \ / -i „ ,sin(7T£).o 7rm 2 ,o , 1 . . 
e(po) (1 + 2ef (^^) 3 (—pT? (^) (3-14) 

The other term {3A 2 B} is a bit more involved. Integrated over P 2 it results in 

e(Po) H^)(~— 2 ) 2 (1 + 2e) - 2(l + 2£)cos(7T£) + (1 + 2e) 2 cos'fre) 1 ^ 

Expanding in powers of e the two trigonometric functions, we get 

^o)3(^)(-^) 2 ^(l + ^ + 0(£ 3 )j (3.16) 

The mass singularity compensation is now made obvious as the first term in the last parenthesis of (3.16), the 
+1, compensates for the real term (3.13), whereas the second term and (3.14) cancel each other, just leaving 
regular contributions. Note that with respect to the previous case, N = 2, expanding the trigonometric 
functions has become mandatory in order to manifest the singularity compensation. 

As we have been able to check that a similar compensation of all leading and subleading mass singularities, 
occur at N = 4 and N = 5 as well, these results strongly suggest that, for all N, such an overall cancellation 
of mass singularities is extremely likely to take place, and not just at N < 2, as was incorrectly stated in [8]. 
However, the control of the ensuing finite parts and their subsequent summation over N, remains somewhat 
puzzling. This is why we will carry out the summation over N a different way. 

First, we take the complicated expressions (2.3) and (2.5) to the much simpler forms 

A(P 2 ,p 2 ) = —m 2 — -y{l — l-yf} , B(P 2 ,p 2 ) = -nm 2 (-^-)^Q(-P 2 ) (3.17) 

£ p p p z 

Of course, with (3.17) as a crude approximation of (2.3) and (2.5), only the leading mass singularities cancel 
out for all N, while leaving a finite series of subleading mass singularities. The latter, though, can be 
resummed into regular functions. Then, bringing the expressions (3.17) to completion, that is to the forms 
(2.3)-(2.5), the full result can be derived rather easily out of those regular functions. This procedure will 
further open two interesting possibilities, dealt with in sections 5 and 6. Also, it applies straightforwardly 
to the (fermionic) cases where HTL self energies are proportional to the second kind Legendre function 
Qo{po,p) rather than Q\. 

We begin with showing that for all N, the stronger singularities, of order e~ N cancel each other. In effect, 
at general N, such a singularity comes from the real part of S^. It is 

((-1) N /N\) 7r£(p ) S^(P 2 ) A N (P 2 ,p 2 ) = ne(p )6(P 2 ) (-m 2 /sp 2 ) N (3.18) 
whereas that very contribution from the imaginary part of can be written 

■^-^-njr.^I^)*^-'!!-^)"-' 
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where the are the binomial coefficients, and where the 0(— P 2 ) distribution has been turned into the 
[— p 2 , 0]-integration range, by taking the condition P 2 « p 2 « T 2 into account. The integral over P 2 
splits into N integrals, the sum of which reads [11] 

<^> (j) (<* | _ ,) ,3.20, 

By putting this expression back into (3.19), one gets the exact compensation of the e~ N singularity of 
equation (3.18), thus completing the proof that the two most singular terms cancel each other at any 
number TV of one-loop self energy insertions along the P-linc. 

It is interesting to note that, for this process at least, the choice of the two independent variables P 2 and 
p 2 is a most convenient one. In particular, singularity compensation is rendered manifest after a unique one 
dimensional integration, the one over P 2 , is performed. Also, it displays the mass (or collinear) singularity 
character of the encountered poles. Had we choosen p and \p\ as another, more customary choice of two 
independent variables, none of these two properties were obtained. 

However a milder singularity (like (3.14)) remains non cancelled. The phenomenon of course starts out at 
N = 3 with the {B 3 } term, and is, thereof, general. For any N = 2k + 1 with k > 1, the imaginary part of 
develops a similar singularity under integration over P 2 

In the imaginary part of S^, but for N- values greater than 2k + 1 at some given value of k > 1, each of 
these terms do appear again, multiplied by some powers of A(P 2 ,p 2 ). Typically, 

e( P o)(-l) k C 2 N k+1 A N - 2k - 1 B 2k+1 , N>2k + 1, k>l (3.22) 

Selecting large enough a value of A, the sum of these terms from Af = 2k + 1 to Af = N, say, reads 

^(po)(-i) fe (-^) 2fe+i E„(-^)^ 2fc - i ^ +i £ 2 dp'p ^-i^g^r 1 ^ 

Carrying out the integration over P 2 and taking advantage of the arithmetical identity (see Appendix) 

2k 1 U ~ 2k ~ 1 • 1 

the expression (3.23) can be given the much simpler form 

r ™? V2k+1 1 

p 2 + 2s ' (2k + 1) e ^M=2k+i 1 £ p 

rt of and letting N tend to infinity, 

the series by singling out the {i3 2fe+1 }-terms we have just analysed. One gets 

^o)(-l) fe i-^) 2k+1 (2feTiJ7 >< Z ^P) > ^) = (l " (-^)) 1 ( 3 - 26 ) 

The occurence of the function Z is quite remarkable as it can be obtained out of the real part of the self 
energy by differentiating it with respect to P 2 , at P 2 = 0. We have, the prime indicating such a derivation 

Z(e,p) = (l-ReE'(0,p 2 )) _1 = (l + m 2 / ep 2 )' 1 (3.27) 
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Summing over AT the imaginary part of and letting A tend to infinity, one can re-arrange the sum of 

,2 



Eventually, summing over k, we get 



oo „ 

TV? 



2fe=JV-l-e((-l) N ) 

dp2 TptWTI E c^+^-l)^ 2 ^ 1 ^- 2 ^ 1 = e(p ) ^>F(p) (3.28) 



N y ' k=l 



where one has 

^Z { s, P )/s = P 2 /m 2 , F(p) = — 2 ^C-r + ') +t g -^^(-Pjf (3.29) 

two regular functions of p in the domain T>. Note that for the series (3.28) to converge and yield the function 
F(p), it is necessary that the condition 

l> \{^-? p>(-^)gT (3.30) 

be fulfilled. At its turn, (3.30) is obviously consistent with the assumed softness of the P-line. The conver- 
gence of the series displayed in (3.25) and leading to (3.27) is less obvious. This is because m 2 /ep 2 cannot be 
expected to lie within the convergence radius of (3.25) when e _1 regularizes a potentially divergent behaviour. 
As should be made clear through sections 5 and 6, though, (3.27) can be provided with the same rigorous 
proof as the one developed at T = 0, relying on the renormalization group equations and the asymptotic 
freedom property of the model [13]. 

To summarize, we have obtained 

Imn fl (Q) = ^-^(^-^(l + 2n B (|)) (Z(e, q/2)/e) F(q/2) (3.31) 

an infrared safe result. Up to the Born term which corresponds to zero self energy insertion, N — 0, the 
result (3.31) should coincide, by construction, to a calculation of Imtl^(Q) in which a bare internal P-line 
has been replaced by the corresponding effective propagator *A R (P). When q is soft, though, this only 
replacement is not able to yield the full order g 2 T 2 correction to the zeroth order calculation, and a similar 
replacement must be envisaged for the second internal line as well. 



Within obvious notations, the imaginary part of the Higgs "polarisation tensor" obtained by inserting N 
and N' self energy corrections along the P and P' lines, see Fig.l, can be given the compact form 

lmn£' N '\Q) = J d D P(l + 2n(p )){lm(A(f' ) (P'))Disc P (A(f ) (P)) + (N - N')} (3.32) 

where ImA^ \P') can be read from equation (3.3) right hand side by replacing P and N by P' and N'. 

Writing the contribution to ImII^ V ' JV \Q) of a generic term of ImS^ \P') (sec (3.3)), one gets a rather 
cumbersome expression which is given in the appendix. The key observation, though, is as follows. Both P 
and P' lines develop singular behaviours on the light cone, but for the P' line, the condition P' 2 = gets 
translated into the condition {2p — q = 0}, and the corresponding light cone singular behaviours into generic 
factors 

e(2p-q)[q(2p-q)]- 1 ^ 2k ' +1 ^ (3.33) 

The point is that integrations over P 2 and then over p, just decouple and get the two mass singular behaviours 
factorised. This amazing simplification is entirely due to the peculiarity of the process under consideration, 
with q = 0, and does not seem to extend beyond. Having integrated over P 2 , and summed over N, one gets 
effectively 

°° j-D— 1 

J2^U^' N \Q) = / l^(l + 2n>))Im(A(f \q-p,p)) (Z(s,p)/e)F(p) (3.34) 

N J P 
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Now the functions Z/e and F of (3.29) are perfectly regular at p = q/2, so that the same pattern as before, 
at N' = 0, applies here again for \P')- Using the P'-adapted form of identity (3.20), that is, explicitly 

4/'^ d-i^i-r'^ <**) 



q/2 [q(2p-q)} l -c V V 1 J 2q e 

it is shown in the appendix that the two most singular terms, of order e~ N compensate each other, while 
leaving a finite series of sub-leading mass singularities. Summing the latter series, the final result reads 

ImIT K M) = ^_^ T (|)^^£(i + 2n s (|))(Z(e, (Z /2)/ £ ) 2 F 2 ( g /2) (3.36) 

It is infrared stable and in line with (3.31) on view of the peculiar {P, P'} internal line symmetry. By 
construction, summing over N and N' from zero to infinity, just replaces both internal lines by their corre- 
sponding effective propagators, *A R (P) and *A R (P'), so that the above result (3.36) has a simple relation 
with the full estimation of the Higgs polarisation tensor H R (Q) imaginary part, by means of the effective 
perturbation theory. Explicitly, 

lm*Il R (q,0)~ lmn| orn (g,0) + 2 x (3.31) + (3.36) (3.37) 

where, within standard notations, the quantity Im*II corresponds to the picturial representation of Fig. 2, 
and where ImII Born is just read from (3.31) or (3.36) by dropping the Z/e and F functions. Once completed 
with the results of the next section, the expressions (3.31) and (3.37) will provide the basis of a future 
numerical analysis [18]. 



4- The full HTL treatment 

Our previous treatment of self energy insertions is based on (3.17), that is on the leading, singular most 
terms in an e expansion of the functions A and B. However, nothing guarantees that the less singular parts 
could be safely ignored. This possibility even appears somewhat paradoxical. Singular pieces, generated 
by this order e _1 -part, usually cancel out in the calculation of transition rates, as we have just seen, and 
are accordingly not expected to be the most relevant ones. Indeed, we are carrying out the calculation a 
different way, and we must now show how our previous results get modified when proper account is taken 
of the subleading ^-expansion terms. For the sake of simplicity this is illustrated by considering zero self 
energy insertion along the P' line. 

As a shorthand notation, we introduce v, the variable 

1 + 2 1+£ e 

v(e) = , or v = l + 2e — > v = 1 + ne + 0(e 2 ) (4.1) 

4 £ 

The second possibility for v corresponds to the case where factors of 2 ±e are dropped, as we have seen in 
the previous section 3, whereas the last expression summarizes any of the two possibilities. The real part of 
the full HTL self energy reads accordingly 

m 2 f P 2 \ 

A(P 2 , P 2 ) = i-^P 2 ) \1 - v\^\ e (©(P 2 ) + cosM9(-P 2 ))j (4.2) 

whereas the imaginary part, the function B(P 2 ,p 2 ) of (2.5), can be written 

B(P 2 ,p 2 ) = -u(e)(-P 2 r^(-P 2 ) , u(s) = v(e) ^ (4.3) 
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Acting upon A N , the terms involving 5^ distributions are left the same as before, whereas the only changes 
come from the ImjE^I-pieces which now read 

2fc=Ar-i-e((-i) lv ) 

Im{S N } = e(po) ]T (-l) k C 2 N k+1 A N - 2k - 1 B 2k+1 (4.4) 

fe=0 

and where the sum over k is extended to k — so as to cover the new singular and regular perturbative 
contributions induced by the new functions A and B. Summing over N , and taking the large N limit, one 
has to control the extra pieces which are generated when e is eventually sent to zero. Once folded with the 
corresponding principal value distribution of (3.3), the building blocks of (4.4) yield 



<«(«))»+' i [<&"t' HM^H ^,) (45) 



where a sign e(po) is understood, as well as the range {1 < iV < co}. Let n(cos(7re)) be that very expression 
between parenthesis in (4.5), and first calculate II(cos(7r£)) at cos(7re) = 1, that is, 11(1). With respect to 
the previous calculation, equations (3.21) to (3.25), the whole change is thus translated into the occurence 
of the new factors u(e), as defined in (4.3), and {w- 7 cos- 5 (7re)} between the parenthesis. Setting cos(7T£) = 1, 
the corrections to (3.21)-(3.25), can be obtained by Taylor expanding 11(1) in powers of e, writing 

fi (i) = E ^ sw (i), fi( ° ) ( 1 ) = 2fcTi (46) 

n=0 

where the value of n(°)(l) is given by (3.24). Unfortunately, when doing so, the coefficients LT^™) result in 
finite series of Stirling numbers of the first kind. While instructive for some respects (next section 6 and 
appendix), the resummation of the ensuing series is practically undoable [23], but for a few (three) Stirling 
numbers. This suggests to follow a different approach. Indeed, a first order inhomogencous differential 
equation in the variable v is readily found for 11(1), that we hereafter define as H(v). Taking the appropriate 
"boundary condition" into account, the solution reads 

n(«) = -i TT (li(l) + C 2k+1 J\xx 2k (l- X f- 2k -^ , 11(1) = ^ (4.7) 

where the value of 11(1), the "boundary condition" is, again, fixed by (3.24). It is clear that the first term 
of (4.7), proportional to 11(1), just reproduces a result identical to (3.31), that is 



(4.8) 



e m 2 ) tg (w(e), 
whereas the second term of (4.7) yields 

^ ( - 1)fc ($y) 2fc+1 \ (5 fe+1 f -i)"-"- 1 ) (+$) N - 2k - 1 ( 4 -9) 

Then, relying on an average value theorem for the integral over x, we know that there exists some number 
c(e) such that (4.9) reads 

^,,H)'(t)" w '^ ^^r'"' ^"s*' <«»> 

To proceed further, it is possible to take advantage of the trick introduced in the appendix, equation (A.l), 
and one finds 
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In order to get (4.11) out of (4.10), series have been summed up, which are convergent within a radius 
specified by the condition 

p A /m 4 > 7T 2 + 2 2 (4.12) 

The limit e = can be taken safely. One has v(0) = 1, and likewise, c(0) = 1 by virtue of (4.10). Adding up 
(4.8) and (4.11), the full result reads 



4) ^ (-r i2 

m z J \p z 



(4.13) 



In the calculation of lmUji(Q), p is fixed at q/2 by the kinematics and so, including the Born term as a 
constant that we do not write for short, one gets eventually 

{0 -}-WI.(0) = »g- (^ ? ) (4.14) 

Now, one must envisage the change in the calculations of (4.5) based on IT(cos(7r£)) instead of 11(1). Indeed, 
one can show that (4.8), (4.11) and (4.14) are left the same, but for terms of order e at most. In the 
calculation of (4.5) based on II (cos(7T£)) instead of 11(1), we get, as a first non trivial correction to 11(1), 

1 \2 I N-2k-l , ,,j . j \ 2 

that is, 

{^}E(-^r +i cr\^)^^ (4.i6) 



N,k 

-2 



The second term of order e is obviously finite. The first one can be calculated and is finite too. One gets 

. 7r 2 , / irm 2 p 2 7TTO 2 \ , 0n , 

2! \(p 2 - Km 2 ) 2 +tt 2 m 4 p 2 J ' 

Then, discriminating between the cases j even or odd (for the expansion of cos^ (tts)), a similar conclusion 
is reached for the subsequent, higher order £-corrections to 11(1). Alternatively, one may proceed a more 
global way, and try to express n(cos(7re)) as a function of 11(1) and cos(7re). A first order inhomogeneous 
differential equation satisfied by II (cos(7re)) can be derived, the solution of which reads (with the appropriate 
"boundary condition", that is 11(1) = H(v)), 



-i / y-i2k+l <>vcos(tts) 

^rr— n( w ) + ^— / dxx 2k (i 



n(cos(7re)) = -^^jrfYTZ-" |n(r)+ / <1, : ^ ] (4.1S) 

Using (4.18), it is straightforward to check that the same result as (4.17) is obtained. 



5- The mixing of singularities 

A well known feature of real time formalisms is the fact that T = and T^O parts come out additive. 
This property of first non trivial order of ordinary perturbation theory does not persist at higher orders which 
mix up zero and non-zero temperature contributions. When the latter are singular, the situation is the one 
of "singularity mixing", an admittedly complicated one [9,12], which, at least at the author knowledge, has 
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not been much investigated until now. In the present situation though, the singularity mixing problem turns 
out to be unexpectedly simplified and this is worth emphasizing. 

At T = 0, the one loop order self energy is well known [9,13] 
Rcnormalized at mass scale /i, one therefore has, writing Hp = a + ib, 

a{p2 ^ 2) = -We p2 i 1 ' j?'" (9(_p2) + e ( p2 ) cos ( 7r£ ))) • b ^ = -if (^ 2 ) 1+e e(p 2 )^^ 

(5.2) 

In the function a(P 2 ,yU 2 ), the first term, given by the +1 in the parenthesis, is the ultraviolet counter-term 
of the T = renormalization procedure. When differentiated with respect to P 2 , though, the same factor 
1/e is to be understood as the dimensional regularization of an infrared (collincar) singularity, as is made 
obvious by simple power counting arguments. 

Now, aT/0 real time formalism can be used to calculate pure T = quantities, just by letting T tend 
to zero. This long known property of real time formalisms [14] has been explicitly verified in [9] at second 
non trivial order of perturbation theory. Thus, in principle, the global one loop self energy correction can 
be inserted in our previous equations, so as to get by the same token, both zero and non zero temperature 
corrections. In practice however, the different analytical properties involved at T — and T^O, have long 
made this property short of any practical purpose. 

Here instead, we will take advantage of the results of previous sections 3 and 4, dealing first with the functions 
A and B of section 3. Then, remarking that D Ul the infinitesimal dilatation operator in the variable u e —o(p), 
ignores both Z(e,p)/e and T — variables, the results are made complete by acting upon them with D u , 

«P (- in (l-^^) B „).g-(«) = «g-(«(l -=£)"') ■ A,-!; M 

Thus, one can re-write those relevant parts of (2.3) and (2.5) the following way 

A{ p2y )=T 2e { ^L f ^_|p2 ^_ |^ r( @ ( p2 )+0( _ p2)cos(7re)) ^ (54) 

and 

B^y) = (^) 2 (-i(-^ 2 ) 1+£ e(-P 2 )^) (5-5) 

A global one loop renormalized (Feynman) self energy can accordingly be defined as 

S F (P 2 ,p 2 , M 2 ) = A{P\p\f) + iB(P\ P 2 ) (5.6) 

with 

A(PW) = a(P 2 ,^) (— ) 2 a(-P 2 ,p 2 ) (5.7) 

P 

B{P\p 2 ) = b(P 2 ) + {^fb{-P 2 ) (5.8) 

Folded with the principal value distributions, over their respective kinematical domains, the functions b(P 2 ) 
and b(—P 2 ) involve the same pole structure, and, up to terms of order (p 2 //i 2 ) £ , the associated residues are 
the same. In these topologies however, these terms do not regularise any singular behaviour and their e = 
limit can be taken with the alluded result (likewise, for the same reason, writing (5.7) and (5.8), the T 2e and 
(T/p) 2e terms of (5.3) and (5.4) have been ignored). 
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As a consequence, the whole series we have just analysed in the pure thermal case, goes through, unchanged, 
in the global treatment of both T = and T ^ parts, but for the two replacements 

a ( l ( ( 27rT )2 



ep 2 12e y \ x p 

^ ^ W _ m 2 (v(e) 1) ™ A + ( M: )2(1 _ I^Kfllll))-! ) 

pz pz £ pz 12 \ p ep z 

Result (4.14), for example, is transformed into 



{C st }- 1 Imn fl (Q) = ( Z(£ ' g/2) j F(g/2) (5.10) 



where we have defined the global functions (adding up the Born term) 



Z-\s, P ) = l + -^-(l + (—f) (5.11) 



12e V v p 

Expressions (5.11) and (5.12) are obviously reminiscent of the elementary first order additivity of the T = 
and T ^ parts, while displaying its full one loop leading order (HTL) realization. Of course, we do not 
expect that this simple pattern should extend beyond the leading HTL approximation. Clearly, the HTL 
peculiarity is at play (see remark (i) below). By construction, the thermal contributions of (5.11) and (5.12) 
are valid for soft internal momenta and are thus enhanced by a factor l/g 2 with respect to the T = ones. 
Some remarks are in order. 

(i) -The striking similarity between the renormalized T — self energy and its T^O, HTL counter part 
is worth emphasizing (compare equations (5.2) with (5.4) and (5.5) ). It is at the origin of a possible global 
treatment, as we have just sketched. Unfortunately, we have been unable to find any simple interpretation 
of this amazing fact, except that it is certainly in line with [21]. 

(ii) -Things can be viewed the other way round. Once it is verified that the couples of functions (a, b) 
and (A, B) are involved in phase space domains over which they develop the same singular structures, then, 
(5.7) and (5.8) obviously indicate that the issues of singularity compensations will be the same in either 
cases. For the topologies considered here, where only mass singularities show up, the above statement is 
nothing but an alternative form of a Niegawa's recent result [15]. 



6. Renormalization constants 

In the past twelve years, efforts have been sometimes devoted to the possibility of defining renormaliza- 
tion constants in a thermal context [16]. When fermionic fields are involved, for example, such a notion could 
be shown to be almost unreliable [16], whereas more room was left in the case of Bose-Einstein statistical 
fields [17]. However, even in the simpler scalar case, some ambiguity was left regarding its definition. 

According to [17], the external leg renormalisation constant should be defined as 

Z(p 2 )=Re(l-S'(0,p 2 ))^ (6.1) 

where the P 2 -derivative indicated by the prime must be identified with a total derivative, kinematics render- 
ing p 2 a function of P 2 , process dependent, though weakly. At T = 0, partial and total derivatives coincide 
because of Lorentz invariance, and the definition of Z reads instead 

Z = (l-RcE'(O)r 1 (6.2) 
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Besides the familiar lack of Lorentz invariance inherent to non zero temperatures, and reflected in the p- 
momentum dependence, definition (6.1) differs significantly from (6.2) by the very place where the real part 
prescription, Re, should be taken. Furthermore, it is worth remarking that, if pertinent, (6.1) should also 
apply at T = for unstable particles [17]. In a few particular cases both definitions agree. For massless 
fields within the dimensional regularization scheme (see (2.3) and (2.5)), for example, and also in the case of 
a real valued self energy function T,(P 2 ,p 2 ), because of the +ie a po advanced or retarded prescription which 
completes the dressed propagator determination. But the one loop self energy is not real valued in general. 

Indeed, the identification (6.1) is obtained by a formal resummation of the terms involving S^ N \P 2 )- 
distributions only, as they appear on the right hand side of (3.3). Now, our analysis has shown that these 
contributions can not be disentangled from those involving principal value distributions, folded with imagi- 
nary parts of self energy insertions. When the latter are non zero, that is, when the self energy is complex 
valued, it was recognized in [17] that an expansion in the number of self energy insertions was necessary in 
order to give a well defined meaning to the calculations. 

The analysis of the previous sections relies on such an expansion and is accordingly able to identify 

Z(p 2 ) = (l-RcS'(0,p 2 ))^ (6.3) 
rather than (6.1), as the correct renormalisation constant, thermal counterpart (see the appendix). 



Appendix 



On section 3 

As the arithmetical identity (3.24) is used extensively throughout the article, we find worthwhile to give 
its demonstration, the more as a proof by induction appears hopeless. Rather we can write 

iV-2fe-l „i 

(3.24) = Cf +1 ( J2 (- 1 ) J '^- 2fe -i) / te* j+2k (Al) 
Now, interchanging the sum with the integral, one obtains the alluded identity as 

C 2 N k+1 / dx x 2k (l-x) N - 2k - 1 =Cf +1 B(2k + l,N -2k) = (A2) 

Jo 2fc + 1 

where B(x,y) is the Euler beta-function. A whole series of similar arithmetical identities can be obtained 
that way. 



Here are given the basic but rather cumbersome expressions leading from (3.32) to (3.36) . From (3.32), 
let us write explicitly 

Im (A<fV)) = P ^gpg - nt^ £ip ' o) S^\ P ' 2 )Re^'(P') (A3) 
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Consider Imll^'^ ' (Q), with N' fixed, and write the contribution of a generic term appearing in ImE^ (P')- 
One gets 



/ 



D _ 1 2k'=N'-l-e«-l) N ') 2 2 



/ 



, 8( P' 2 )P , _ I -P' 2 |£NjV , -2fc , -l /__-XCp2w_^L.\iV P 

(_p/2)l-(2fc'+l)e V 1 1 4p/2 1 - 1 | ""^ A £p 2> (p2)JV+l 

2fe=JV-l-e((-l) N ) 2 . 2 \ iV-2fe-l 

(- 1 ) fe (-^) 2fc+1 ((- p2 ) 1+£ ) 2 ' £+le (- p2 )} ( M ) 

In the light cone neighbourhood, the sum over P 2 develops mass singular behaviours , whereas the remainder 
of the integrand is a well behaved function. The singular part of (A4) is accordingly given by 



/ 



^(i +2 n s ( P ))^;+ i (-i) fe ' e ( 9 -p)(-^r+ i (-^r- 2fe '- 1 

N'-2k'-l 



-) N 

2 > 



( 



-i) k (-0;) 2k+1 ((-p 2 ) 1+e ) 2k+1 e(-P 2 )} (A5) 



where we have used the relation p = jf, particular to both the physical process under consideration with 
Q = (q, 0), and our notations, specified on Fig.l. Integrating now over P 2 , the singular most terms of order 
e~ N cancel out thanks to (3.20), leaving a finite series of sub- leading mass singularities. The latter can be 
summed over N, using (3.24 ), so as to yield eventually 



D _ 1 2k' = N'~l-e((-l) N ) 2 

£(^) = / V^ (1 + 2nB(rt) £ ^\-lfe(q- P )(-^^ 

N J P fc ,_ Q P 



The integral over p factorizes the total solid angle of the model. Then, the integration over \p\ ranges from 
p = q/2 to p = 00 because of the <d(2p — q) distribution, but the sign distribution e(q — p) splits the integral 
into two parts. The second part, say, with an integration range from p = q to p = 00, will be discarded 
because, then, one is leaving the domain of softness of the P, P' lines, that is, the realm of relevance of the 
effective perturbation theory. Even though, it is easy to check that this neglected integral is perfectly regular. 
The first part, with an integration range from p = q/2 to p = q, exhibits light cone singular behaviours at 
the p = q/2 lower bound, through the factors (3.33). Now, by using (3.35), it is straightforward to show that 
the singular most contribution, of order e~ N x (q 2 /4m 2 ) x F(q/2) coming from the term k' = in (A6), 
cancels out against the contribution coming from the ReE^ (P') term of (A3) exactly like for the P line 

-« { -^<Po)5 {N '\P' 2 )Rexf\p<) => -Mp'o)S(P' 2 ) (-^Y (AT) 

In (A6), one is left with the sum of {P 2fe +1 }-terms, at k' > 1, which entail the factor {(q(2p — q))~ 1+( - 2k 
singular at the lower boundary of the p-integration, whereas the remainder of the integrand, that is, basically, 
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the functions Z(e,p)/e and F{jp), is perfectly regular. Then it is easy to use the arithmetical identity (3.24) 
in order to show that the same pattern as for the P-line applies here again, with a finite series of sub-leading 
mass singularity outcome. Summing the latter over TV', equation (3.36) results in the limit N' = oo. 



On section 6 

In section 3, equation (3.26), we have pointed out the occurence of a global factorizing expression Z(e,p). 
Together with the overall factor l/s, the function Z contributes the global simple function p 2 /m 2 (as seen 
from (3.29)), which is recovered also throughout the analysis performed insection 4. In this latter section 
though, the function Z is not itself immediately apparent, thereby obscuring somehow its interpretation as 
a renormalization constant. It is indeed an artefact of the approach followed in section 4. This is realized 
by getting back to the e-expansion (4.6), for which we have (N — 2k — 1 > n) 



»->(!) = csr "if ' ( ~ 1)J, ";a"r" +1) c ^- < A9 > 

Now, by an iterated use of of the arithmetical identity (3.24), and also of the identity [11] 

N-2k-l 

]T j C j N _ 2k _ x = 0, 1 <*<iV-2fc-l (A10) 

3=0 

it is easy to show that any of the n(") (1) is a polynomial of degree (n — 1) in the variable (2k + 1). Explicitly 

n-l 

ft n )(l)=Y J (-l) m S™_ 1 (2k + l) m , n>l (All) 

where the coefficients {S^ l _ 1 } arc integers, namely, Stirling numbers. This can be realized by quoting the 
relation, due to (A. 10), 

jj(«)(i) = _[(„_ l) + (2k + 1)] n("" 1 )(l), n>\ (A12) 

and this yields 

n (n) (l) = (-!)"{ II (m+[2k + l})\ (^—) (A.13) 

where the generating relation for the Stirling numbers of the first kind is manifest [23] . 

The point is that at any given order e n in the expansion (4.6), the coefficient IT™)(1) is independent of N. 
On view of (4.5), we conclude that the summation over N can be carried out, and that the ensuing expression 
{Z(e,p)/e} keeps entering the result, for all n and all k, as an overall multiplicative factor. 
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7. Conclusion 

On the basis of a scalar model of field theory, renormalizable and asymptotically free, the light cone 
singular behaviour of the corresponding thermal effective perturbation theory has been investigated for some 
given process. 

Using the HTL form of the self energy near the light cone where mass singularities develop, one can conclude 
that, at any number of self energy insertions, a full compensation of singularities is obtained. This conclu- 
sion brings a correction to the erroneous statement which appeared in [8], and was due to an insufficient 
dimensional expansion of the self energy functions. In this respect, it may be recalled that after a first 
series of calculations performed at first and second non trivial orders of thermal perturbation theory [3,9,19], 
it has been observed (and then conjectured) that singularity cancellation took place within each topology 
separately [20]. We can remark that our analysis supports and enforces this conjecture. 

Once the detailed balance compensation of mass singularities is achieved, the result for the process under 
consideration, can be obtained by summing over {N, N'} the regular contributions which are left, as this 
summation is mandatory in order to get the full leading order correction 0(g 2 T 2 ) to the Born approximation. 

Throughout the present analysis, however, the summation has been carried out a somewhat different way. 
Although the full HTL self energy is recognized a perfectly regular function of P 2 in any light cone neigh- 
bourhood, we have stressed that it entails a potentially singular most part, dealt with in section 3, and which 
turns out to yield the basic result. Remarkably enough, at the level of approximation where we have been 
working, this basic result enjoys two interesting properties. First, it is immediately related to the cases in 
which the HTL self energies are given by the (second kind) Legendre function Qo (fcrmionic fields) and is 
therefore relevant to the study of thermal QED and QCD, the infrared problem included [22]. Second, it is 
simply enough related to the case under consideration, where one has to cope with a HTL self energy on the 
order of the Legendre function Qi rather than Q . 

Furthermore, proceeding that way, a striking analogy can be put forth between the T = rcnormalized 
self energy and its HTL counterpart. This is certainly a remarkable and intriguing aspect which, at least 
in our opinion, can only be read in line with the ideas developed in [21]. In particular, the issued mass 
singularity structures (poles and residues), turn out to be identical, and a global treatment of both zero and 
non zero temperature corrections is made possible, resulting in a simple resolution of the otherwise intricate 
singularity mixing problem. Likewise, a new, unambiguous identification of a rcnormalization constant in 
the thermal context is obtained, and turns out to be at variance with some previous attempts. 

The function F(p), and a fortiori, its generalization F(p), is new and interesting in this context, as it could 
not be obtained in previous related works [12], [18], where the effective perturbation theory was used right 
from the onset. Remarking the similarity of F with a first order brehmsstrahlung function, one recovers a 
physical picture of the HTL's which is closely related to their interpretation in terms of forward Thompson 
scattering amplitudes [21]. 

Eventually, we can stress that our guiding strategy has been to express the singular structure of an effective 
quantum field theory at temperature T, in terms of the same bare temperature quantum field theory. Such 
a strategy may be thought of as arduous. On the other hand, though, the complexity of the effective 
perturbation theory, the mixing of topologies it envolves, renders hazardous the control of such fine tuning 
mechanisms as singularity compensations, if not the very meaning of the ensuing finite results. The same 
approach is extended to the case of gauge field theories for timelike and lightlike external momenta [22]. 
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